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Spin-orbit (SO) interaction critically influences electron spin dynamics and spin transport in bulk 
semiconductors and semiconductor microstructures. This interaction couples electron spin to dc and 
ac electric fields. Spin coupling to ac electric fields allows efficient spin manipulating by the electric 
component of electromagnetic field through the electric dipole spin resonance (EDSR) mechanism. 
Usually, it is much more efficient than the magnetic manipulation due to a larger coupling constant 
and the easier access to spins at a nanometer scale. The dependence of the EDSR intensity on 
the magnetic field direction allows measuring the relative strengths of the competing SO coupling 
mechanisms in quantum wells. Spin coupling to an in-plane electric field is much stronger than to 
a perpendicular field. Because electron bands in microstructures are spin split by SO interaction, 
electron spin is not conserved and spin transport in them is controlled by a number of competing 
parameters, hence, it is rather nontrivial. The relation between spin transport, spin currents, and 
spin populations is critically discussed. Importance of transients and sharp gradients for generating 
spin magnetization by electric fields and for ballistic spin transport is clarified. 

PACS numbers; 71.70.Ej,72.25.Dc,72.25.Hg,78.67.Dc 
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INTRODUCTION 

Manipulating electron spins at a given location and 
transporting electron spins between different locations 
belong to the central problems of semiconductor spin- 
tronics and are of critical importance for quantum com- 
puting and information processing Q, |^ |j, Q . Among 
the different concepts of spin injection and spin manip- 
ulation that are discussed in the current literature, the 
approaches based on SO coupling acquire growing atten- 
tion. 

Spin orbit interaction couples electron spins to the elec- 
trical component E (t) of electromagnetic field. The 
Hamiltonian in many cases can be presented in a form 



H{r , k, a) = Horbi k, (t) + Hz{r , cr ; 



(1) 



where r , fc , and <r are coordinates, momenta, and 
Pauli matrices, respectively. Each term of the Hamil- 
tonian including the operators of spatial quantities ( r 
and/or k ) and Pauli matrices produces SO interaction. 
In Eq. (^1, the first term Horb{ k , a) symbolizes orbital 
mechanisms of SO coupling dependent on the electron 
momentum, while the second term Hz ( r , cr ) symbol- 
izes the Zeeman energy for a system in an inhomogeneous 
magnetic field or with a spatially dependent -factor. 
In the bulk, orbital mechanisms of SO coupling usually 
lay the major role |^ Q and result in a strong EDSR 
Ig] dominating over the electron paramagnetic reso- 
nance (EPR), the magnetic excitation of spin transitions. 
However, the spatial dependence of the Zeeman interac- 
tion may also play role 0, 0, 0| . Remarkably, electri- 
cal operation of electron spins in Al ^ Ga i^x As quantum 
wells has been achieved first through the Zeeman mech- 
anism. Kato et al. took advantage of the anoma- 
lously small g -factor of the bulk GaAs, g « —0.4 , that 
allowed them to achieve both the strong anisotropy of 



the g -tensor and its spatial dependence across an inho- 
mogeneous quantum well, g — g{z) . Meantime, orbital 
mechanisms of SO coupling are also greatly enhanced in 
quantum wells as compared with the bulk material be- 
cause of the lowering the symmetry that results in devel- 
oping new SO terms in the Hamiltonian. It is shown in 
Sec. II that these mechanisms should result in a strong 
EDSR, especially in materials with large g -factors typi- 
cal of narrow-gap semiconductors |l2Lll3|| . Extraordinary 
efficiency of an in-plane field E (t) is one of the basic 
conclusions of Sec. II. 

After a number of exciting proposals about using SO 
coupling for electrical spin injection, it has been already 
achieved experimentally |14lll5l |. This success pushes for- 
ward the problem of the propagation of inhomogeneous 
spin populations, say wave packets [l^ . in media with a 
SO split spectrum [l7|. It sounds tempting to apply to 
this problem a notion of spin currents defined similarly 
to electric currents. E.g., electrical spin injection and 
spin polarization produced by SO mechanisms and spin 
transport in such systems have been discussed recently 
in terms of spin currents driven by ac 0, 0| and dc 
0i|23| electric fields. However, this appealing approach 
meets serious problems. First, in media with SO inter- 
action spin currents are not conserved, hence, there is 
no rigorous definition of them. Second, there is no ex- 
perimental procedure for measuring them. Third, these 
currents do not vanish even in thermal equilibrium when 
there is no physical spin transport. Fourth, as distinct 
from electric charges and electric currents that appear in 
Maxwellian equations explicitly, only the spin magneti- 
zation can be included into the total magnetization M 
in a straightforward way; the place of spin currents in 
macroscopic physics is still to be unveiled. In addition, 
the properties of these currents found in Refs.Eil andl2ll 
were quite unexpected and inspired an active discussion. 
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Consensus about this problem has not been achieved yet. 
There are interesting attempts to find the Umits within 
which the concept of spin currents, however not rigorous, 
can be appHed. The diff'erent option is to concentrate on 
the magnetization, the quantity that can be both rigor- 
ously defined theoretically and accessible to experimental 
control. In Sec. Ill, some of the related problems are dis- 
cussed. 

SPIN DYNAMICS IN QUANTUM WELLS 

Quantum wells are usually noncentrosymmetric, and 
the two-fold spin degeneracy of the electron energy spec- 
trum is lifted by two SO terms 

Hr = aR{axky ~ cTykx), Hd = aoia^kx - (Jyky). (2) 

Hfi is due to the structure induced asymmetry (SIA) 
and is known as Rashba term, Hjj is due to the bulk 
induced asymmetry (BIA) and is known as Dresselhaus 
term, and ay are Pauli matrices, and k^ and ky are 
the components of the momentum in a magnetic field B 
[2^. Because the basic results are similar for Hr and 
Hd , only equations for Hr will be presented in what 
follows. Intensity of EDSR is controlled by the interac- 
tion e(rso • E (t)) , where r so is a properly defined 
SO contribution to the coordinate operator. When B is 
strong enough, SO coupling can be treated as a perturba- 



where if and ip are azimuths of B and E , respectively. 
With aR w 10"^ eV cm, m « 0.05mo , and B « 1 T, 
we arrive to ij^ as large as « 10~^ cm. Then the Rabi 
frequency Qr = eEl/Ti , that is the basic figure of merit 
of the spin operation efficiency, equals il/? w lO^*' s^^ in 
a field as small as only about E « 0.6 V/cm. With such 
a high efficiency, electrical spin operation seems promis- 
ing even for Si quantum wells despite of the small aR 
values typical of them, aR « 10~^^ eV cm [j^j . 

Equations © and (@J were derived for a strong field 
B and a perfect quantum well. For a moderate field 
B , spin-flip frequency depends on the Landau quantum 
number, and a random potential results in an inhomoge- 
neous broadening. Under general conditions, the EDSR 
band consists of two components with the widths about 
Tp"^ and , Tp and being the momentum and spin 
relaxation times, respectively |2J|. However, when v so , 
the SO part of the velocity operator, does not depend 
on the momentum k , the total oscillator strength of the 



tion, and, for an electron confined in a parabolic quantum 
well, the matrix elements of r so between spin-up and 
spin-down states can be found. 

In a perpendicular electric field^^ E || z , the matrix 

element of a spin-flip transition is |l2| 

/t r ||\_ uj cU> s {uJ c - (^s) sin 20 

{ I |2so| I) - ^2^2 cos2 Q _ ^2(^2 _^ ^2 _ ^2) " 

(3) 

Here loq is the parabolic conflncmcnt frequency, Uc — 
uJcifi) and ujs are the cyclotron and spin frequencies, re- 
spectively, and 9 is the magnetic-field polar angle. When 
both SIA and BIA contribute to (t l^sol i) , it acquires 
an azimuth dependence, and the angular dependence of 
(T l^sol i) can be used for measuring the ratio aR/ao ■ 
With aR ^ 10^^ eV cm (as typical of InAs quantum 
wells), (t |zso| i) is usually considerably larger than the 
Compton length, "Xc — fi/rriQC « 4 x 10"^ cm, that 
plays the role of a characteristic length for EPR. There- 
fore, EDSR is stronger than EPR. 

Nevertheless, the matrix element (t l^sol i) includes 
a factor lUs/ujq that is usually small and reduces the 
EDSR intensity, loq in its denominator indicates that 
the deviation of the system from the strict 2D limit is the 
critical condition for EDSR, and it vanishes for 9 = 0, 
i.e., EDSR can be observed only in a tilted field B. 
Measurements with an in-plane electric field, E ± z , 
allow one to get rid of all these problems. A matrix 
element similar to (T |^so| i) equals ^3 



(4) 

I 

spin-flip transition is contained in the narrow band. Be- 
cause i} so does not depend on k for both Hr and 
Hd , a strong dynamic narrowing of the EDSR band is 
expected. 

In conclusion, EDSR seems to be a highly promising 
tool for the electron-spin manipulation in quantum wells. 



SPIN TRANSPORT IN MEDIA WITH 
SPIN-ORBIT COUPLING: SPIN FLUX AND 
SPIN CURRENTS 

Spin currents are widely discussed and are a subject 
of excitement and controversy. In what follows, my cur- 
rent understanding of the problem is substantiated and 
summarized. 

The notion of the spin flux, or, what is the same, of a 
transport spin current, has an appealing physical mean- 
ing as the propagation of an inhomogeneous spin po- 



OIR 



[{ijJc COS 9 -f LOs) cos((/3 — V') + *(t^c + cos 9) sin((^ — ?/;)], 
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larization. However, a rigorous definition of it is lack- 
ing. Usually, the definition of the flux density J ^ ( r ) 
of a physical quantity A originates from the continuity 
equation for A . The equation j ^( r ) = r ) d ( r ) , 
where p/i ( r ) is the density of A and v {r ) is the ve- 
locity of the flow, holds for the mass flux and electric 
current due to the mass and charge conservation. Be- 
cause the mass m of a particle and its electric charge 
e do not depend on the momentum k , this definition 
is tantamount to the integration of the mass and charge 
currents, mv (k) and ev (k) , over the k -space with 
an appropriate distribution function. However, simple 
relations of this sort are not universal. 

E.g., because the energy conservation law includes the 
work performed by the pressure P( r ) , the energy flow 
density in an ideal fluid equals j ^{r) = pw{r)v{r) 
and includes the enthalpy density Pw{r) — p^lr) + 
P{r ) rather than the energy density Pe(r) 25]. 

In media with SO interaction, spin is not conserved. 
Spin dynamics in the effective magnetic field B cff ( ) 
and, hence, spin nonconservation, is central for the Datta 
and Das spin-transistor concept '2&\. Therefore, it is well 
understood that the definition of the Hermitian operator 
of the spin- current density 

Jrj ( fe ) - \g,V, {k)+Vj{k )(T,]/2, (5) 

that is similar to the mass and charge currents and is 
widely used, lacks the p rop er j ustification for systems 
with SO coupling 123, 123, l22l ; here i,j are Cartesian 
coordinates. To the best of my knowledge, no experi- 
ments for measuring the quantities defined by Eq. Q 
have been proposed. 

Another problem with the spin current is related to 
its behavior with respect to the time inversion, i.e., the 
t —t transformation. The mass, charge, and energy 
densities are even (or real in Wigner's terminology) with 
respect to the i -inversion, while v is odd (imaginary). 
Therefore, the corresponding currents are t -odd and, 
after the integration over k , vanish in thermodynamic 
equilibrium. On the contrary, Pauli matrices cr are t - 
odd, hence, spin currents are t -even and their mean val- 
ues do not necessarily vanish in equilibrium |28l |. Ex- 
istence of such currents in noncentrosymmetric systems 
is in agreement with general principles of statistical me- 
chanics 1^3 ■ 

These observations indicates existence of a gap be- 
tween the physical notion of a spin flux ( = transport 
spin current) and the formal definition of spin currents 
through Eq. JSJ. 

Momentum Current, Momentum Flux, and 
Knudsen flux 

In this context, it is instructive to begin with con- 
sidering first a classical quantity, the momentum cur- 



rent Uii = T,ipi{l)vi{l) that also is even with respect 
to t -inversion. Here I numerates particles inside a 
unit volume, and Pi{l) and Vi{l) designate their mo- 
menta and velocities. For an equilibrium gas Hjj = P , 
the gas-kinetic pressure that is a thermodynamic rather 
than transport property. A similar conclusion follows 
from macroscopic arguments j25j . The mass-flux density 
j mi'"') of a, streaming fluid equals to its momentum 
density p{r) v (r ) , p{r) being the fluid velocity. Mo- 
mentum is conserved in an ideal fluid, and the continuity 
equation d{pvi)/dt = —dUij/dxj results in the equa- 
tion Hy = PSij + pviVj for the momentum-flux density 
tensor H^^ ( r ) ; in equilibrium, its diagonal components 
reduce to the momentum current P found above. The 
nondiagonal part [i ^ j) of the momentum-flux density 
p{ r )vi{ r )vj{ r ) , describing the momentum transport, 
can be found from H^ directly; it vanishes in equilib- 
rium. However, splitting Hij(r,i) into the pressure P 
and the momentum flux requires solving the set of the 
fluid-mechanics equations and cannot be performed in a 
general form. 

Remarkably, for a Knudsen gas flow across a nar- 
row opening, the flux can be easily related to P . In- 
deed, the momentum-flux density equals to the integral 
of the momentum-current density H^^ ( k ) over a single 
hemisphere, hence, the Knudsen momentum flux equals 
Iiii/2 = P/2 , a half of the momentum current. For the 
Knudsen flux, a large gradient of the concentration is 
critical. 

Spin Currents in Thermodynamic Equilibrium 

The similarity between Jij ( fe ) and H^- ( k ) discussed 
above suggests existence of nonvanishing equilibrium spin 
currents in noncentrosymmetric systems; they will be 
termed as background spin currents in what follows |28l | . 
For the Hamiltonian Hn , a calculation renders 

Jxx = Jyy = 0, = ~Jy, ^Jr^ m^ayiTTfC', (6) 

whenever the electrochemical potential /i > . Spin cur- 
rents carried by separate electrons are linear in . The 
current Jr is proportional to a^j because of the par- 
tial cancelation of the contributions from two spectrum 
branches, e\{k) = h^k^/2m + XaRk , A = ±1 . When 
p ^ Ea , where = ma\/TT? is a characteristic SO 
energy, the difference is relatively small. Nevertheless, 
it persists under the conditions of thermal equilibrium 
when there is no spin transport. The nonvanishing equi- 
librium spin current should be taken as a warning indi- 
cating problems inherent in the spin current concept. 

The real antisymmetric pseudotensor J7ij is equiva- 
lent to a real vector P || z . This relation implies the 
existence of some connection between the background 
spin currents and the SO contribution to the polariza- 
tion of the 2D electron gas by the electric field E ±^ 
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producing the SIA. If to define the polarization P as 
P = —d£ii/dE± , where £u is the SO energy of the 2D 
electron gas, then 



P = (Am/h) Jfl daR/dE_L. 



(7) 



Therefore, spin current Jr is directly related to the SO 
contribution P to the electric dipole moment. Here E± 
includes not only the field inside the confinement layer, 
but also the change in the band offsets that strongly in- 
fluences a_R [13 ■ 

Background spin currents are highly sensitive to the 
symmetry of the system. Numerical calculations suggest 
that an can show an essential k dependence. Accepting 
Q!_R(fc^) = ao + ai{k'^) with Q:i(fc^) <C oq , one comes to 



1 + 2fih'^ /maQ 



(8) 

where and k_ are Fermi momenta for the upper 
and lower spectrum branches, respectively. When writ- 
ing Eq. ((SJ), the contribution proportional to [ai(fc?.) -|- 
Q!i(fc^)]/2 has been subtracted; it can be found by dif- 
ferentiating Eq. (jSJ over an . For a small ao , when 
mal/h? <^ fi = h'^kp/2m, Eq. © simplifies 

J,y = (mfcj.|ao|/47r?i3)[a,(fc2 ) _ ai(fc2 )]. (9) 

As distinct from Eq. (O, this expression is quadratic 
rather than cubic in SO coupling. This is a result of 
the symmetry breaking. Indeed, the Hamiltonian Hr 
with aji =const possesses a hidden symmetry that man- 
ifests itself in the fact that the square of the veloc- 
ity operator v = h'^dHn/d k is related to as 
— 2Hji/m + 2a^/h'^ . As a result, its eigenvalues 

— 2e/m + 2a|j,/?i^ depend only on the energy e and 
do not depend on the branch index A . A similar relation 
holds for the semiclassical velocity, v^^ — 2e/m + aj^/h^ . 
Both equations result in A -independent Fermi velocities. 

The Hamiltonian Hr is typical of (0,0,1) 2D layers of 
cubic crystals. For a 2D layer of an uniaxial crystal like 
CdS with an in-plane pyroelectric axis c || y , a new SO 
term 



(10) 



develops in the Hamiltonian. When a^. <C an , spin 
current produced by this term equals |84| 



Jzx = -kpaz/2Trh. 



(11) 



Therefore, in this case the background spin current J^zx 
is linear in the SO coupling constant . Background 
spin currents are confined inside a sample. In a texture 
consisting of zinc blende and wurtzite crystallites, the 
current J^zx is confined inside the latter ones. In equilib- 
rium, the boundaries adjust to ensure this confinement. 



It is seen from Eqs. (O ~ Hll|l that, as distinct from 
Hij , the nondiagonal components of J^ij do not van- 
ish even in equilibrium, and the lower is the symmetry 
the larger the background spin currents are. In a differ- 
ent way, equilibrium spin currents were found by Pareek 
|35l |. Similar problems arise in the theory of magnetic 
materials in inhomogeneous magnetic fields |36Ll37l| what 
is quite natural because the Hamiltonian Hz {r , a) of 
Eq. represents one of SO coupling mechanisms. 



Response to a Time-Dependent Electric Field 

The paper by Sinova et al. |2l| has attracted recently 
an active interest to the spin current J^zx driven by a 
homogeneous electric field E = Eyf . The operators 
entering the corresponding Kubo formula are 



hky/m + aRa^/h, Jzx ^ fi-kxC^z/m. 



(12) 



Because in the basis of the eigenspinors of the Hamil- 
tonian Hn (the chiral basis) the operator jzx has 
only interbranch matrix elements, the spin conductivity 
'^zxy = JzxjEy of a perfect crystal originates from the 
interbranch electronic transitions inside the ring in the 
fc -space restricted by the Fermi momenta k± , Fig. 1. 
The real part of it. 



Kxyi^^) 



eaR I 
2nnm Xc 



k^ dk 



{2aRk/h)^ -LU^ 



(13) 



in the w — > limit results in Yi'^^y(uj = 0) = e/Airh in 
agreement with Ref. EH. Remarkably, T^'^^y{u!) is related 
to the SO contribution, £30(0^) , to the real part of the 
dielectric function e{uj) by a simple equation 



Kxvi^) 



(h^uj'^/STTemaR) 



eso(w)- 



(14) 



Hence, for the Hamiltonian Hr the spectrum- 
specific frequency dependence cancels from the ratio 
'^'zxyi^) / '^'soi^) ■ Eqs. O and ((HI) indicate existence 
of an algebraic connection between spin currents and 
the SO contribution to dielectric polarization. These 
quantities should be related in some way because in 
Maxwellian electrodynamics the response of any homo- 
geneous medium to a homogeneous electric field is de- 
scribed by the dielectric function. 

Deriving E^^(w = 0) in this way from the Kubo for- 
mula of Eq. 1)13(1 is equivalent to the effect of the SO 
term H!,^ with az = —eEy/3kp . The same equation 
can be derived by using the perturbation theory in the 
operator V — —eEyy — ieEyd/dky in the 2x2 spinor 
space. It results in a new equilibrium state in the basis 
of new eigenspinors. From this standpoint, "E'^r^y^uj = 0) 
of Eq. 1(14(1 is a thermodynamic rather than a kinetic pa- 
rameter and is similar to the pressure P of Sec. HI. A. 
Such static spin currents are dissipationless because they 
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do not transport spins. This conclusion suggests that 
^'zxy = 0) calculated by using only interbranch matrix 
elements should be nearly insensitive to impurity scatter- 
ing. Of course, one has to bear in mind that the similarity 
between the two above swtems is not complete. In the 
system by Sinova et al. |2l| there is an external homoge- 
neous electric field Ey , while in a conducting pyroelectric 
the macroscopic field is screened and Ey should be un- 
derstood as an equivalent microscopic field changing the 
electron energy spectrum. 




Momentum 



Fig. 1 Electronic transitions contributing to spin currents 
are shown by arrows inside the gray strips. ^ is the electro- 
chemical potential, k+ and fc_ are the Fermi momenta for 
the upper and lower spectrum branches, respectively; T = . 
Minima of the spectrum e{k) are achieved at the points 
±ka , where ka ~ moiB./h? . 

The problem is that the singular intrabranch contri- 
butions to q = responses (whose explicit form is 
sometimes elusive) nevertheless exist and manifest itself 
when the momentum conservation is broken by impuri- 
ties. E.g., the low- frequency contribution of the Hamilto- 
nian Hj^ to the imaginary part of the dielectric function, 

e^Jc.) = {^uoll2uj)5{u:), (15) 

comes through the rcnormalization oj^ = ^^[1 — 
{maji/'h?') /2TTn] of the "plasma frequency" Wp = 
Ane^n/m, n being the 2D electron concentration [28j |. 
This contribution manifests itself in magnetic resonances 

m 

Generating Spin Fluxes and Spin Polarization 

When a system is out of equilibrium, the states with 
non-zero spin currents have potentialities for generat- 
ing spin fluxes. A nonequilibrium state prepared by an 
electric-field pulse E {r ,t) inside a narrow spot can 
serve as an example. If the boundaries of the spot are 
narrow as compared with the electron mean free path 



vpTp , vp = hkp/m being the Fermi velocity, the density 
of the Knudsen-type spin flux across the boundary can be 
as high as about T,'^^y{uj = 0)E/2 , cf. Sec. III. A. Appar- 
ently, spin fluxes observed by Stevens et al. in their 
optical experiment were of this sort. Strong nonequilib- 
rium near the edge of the spot is critical for generat- 
ing spin fluxes across it. Knudsen spin fluxes across an 
abrupt boundary of zinc blende and wurtzite crystallites 
(Sec. III.B), when the pyroelectric field changes faster 
than the boundary equilibrium establishes, can serve as 
an example. Similar spin fluxes have been discussed by 
Mal'shukov et al. from a different standpoint . 

Electric field E is t -even while magnetization M is 
t -odd. Therefore, for an electric field to produce a mag- 
netization (a magneto-electric effect), t -symmetry should 
be broken by a pr oper magnetic structure [3i||, electron 
scattering |40|, |41| , finite frequency, circular polarization, 
etc. E.g., a dissipationless spin current Jzx driven by the 
field E = EyY is not accompanied by any magnetiza- 
tion, however, a magnetization M \\ x proportional to 
Tp does develop in this geometr y 143 . In microstructures 
with a Sharvin-type resistance |43l |. the resistance quan- 
tum h/e'^ , h — 2Trh , plays the role of the time-symmetry 
breaking parameter. E.g., the resonant-tunneling-diode 
spin-filters proposed by Voskoboynikov et al. |44| are 
controlled b y j Vi while their modification proposed by 
Koga et al. |43 might be controlled by h/e"^ . 

When Tp <C Ti/aiikp , the magnetization is long liv- 
ing due to the Dyakonov - Perel process ^] , EDSR line 
is dynamically narrowed |24j , quasiequilibrium in the or- 
bital degrees of freedom is established 0, 0| , and prop- 
agation of the magnetization M {r ,t) is controlled by 
diffusion, drift, and appropriate SO corrections. With in- 
creasing Tp , different scales related to SO coupling come 
into the game. Spin fluxes not accompanied by magneti- 
zation have short decay times about Tp . 



Response to a Space- and Time-Dependent Electric 

Field 

Therefore, transients and spatially inhomogencous spin 
populations are critical for spintronics, and spin polariza- 
tions are the principal measurables, i.e., quantities acces- 
sible to experimental control. 

Spin transport has two major aspects: generation of 
nonequilibrium spins and their relaxation. There is an 
extensive literature on the mechanisms of spin relaxation. 
In this section, I find the spatial and temporal scales 
related to spin production in media with SO coupling. It 
can be done by calculating the response of free electrons 
to a field E exp{i[(<7 ■ r) — cot]} . The static M^- 
magnetization driven by a field E = Eyy has been 
recently observed 0>S|- 

The in-phase spin polarization in the frequency range 
of the interbranch transitions comes from the S -function 
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parts of Green functions. For q = qk and E 
the mean value of the Pauh matrix Ux reads 



Eyf, 



(Ta:{q,Uj) = 



cos 



(1 + gCOS(/)/2fca)2 



(16) 



Here ka — man/fi^ is the shift of the minima of the e(fc) 
curves of Fig. 1 from the origin. The integration over the 
angle (j) , the azymuth of fe , in the response function 
I{q,Lu) is performed over the interval G $(<7,a;) se- 
lected by equations 



4k + q/2)>fi, £_(fc-q/2) </i, 
£+(fc + q/2) -e-{k- q/2) = hto. 



(17) 



The family of integration paths over </> , parameterized 
by the transition frequency uj , fills the integration area 
in the fc -plane. For q,ka kp , this area is close to 
a ring restricted by two circles with the radii kp ± ka , 
whose centers are displaced from fc =0 by ztq/2. 
The reactive part of spin polarization, phase shifted by 
tt/2 with respected to the field, has the same order of 
magnitude but its analytical form is more involved. 

The position of the band center, iD — 2vpka , and 
the half- width of the band, lji/2 = vp{q + 2ka/kp) , 
determine two characteristic scales in the frequency do- 
main. For small q , q ka , the band is narrow, 
'^1/2 ^ however, — > when q —>■ 2ka ■ The 
topology of the integration paths in the fc -plane changes 
at qcr ~ '^ka/kp . For q < q^r , the shape of the band 
evolves smoothly. Fig. 2. ki q = q^r , the band ac- 
quires a sharp peak, and for q > q^r it becomes wide 
and strongly asymmetric. When q — > 2ka , two circles in 
the fe -plane touch to one another, and the low-frequency 
wing of I{q,oj) moves to w = and diverges. There- 
fore, 2ka/kp and ka are two characteristic SO scales 
in the momentum domain. The divergence of /(g, lo) 
for q — > 2ka underscores the importance of the elec- 
tric field gradients for the production of nonequilibrium 
spins. For q > 2ka , the coUisionless regime is no more 
self-consistent for an infinite area. 

Another contribution to (T^(<7,lj) comes from intra- 
band transitions. When q kp , the in- phase magneti- 
zation does not vanish for lj < vpq , i.e., 



ffxi g , w) oc 0{vpq - Lo), 



(18) 



where Q{x) is a Heaviside function. Similar equations 
exist for spin currents Jij , while their specific form de- 
pends on the indeces i,j . Therefore, the intrabranch 
contributions to the magnetization and spin currents ex- 
ist, should compete with the interbranch contributions, 
and are singular in the q,uj ^ limit. This behavior is 
reminescent of the Schliemann and Loss discussion of the 
dependence of spin currents on the ratio of lo and l/rp 



p9| . It is obvious from Eq. (|18|l why the intrabranch 
contribution to spin magnetization disappears from the 
q = responses in the clean limit. 




Dimensionless frequency, Q 

Fig. 2 Response function I{Q, O.) for the dynamic spin 
magnetization vs dimensionless frequency Q = m(Lo — 
u)) /2TikFka for four values of the dimensionless momentum 
Q = q/2k^ ■ kc^/kp = 0.1 , Qcr « 0.1 . 

Eq. H18|l has been derived for the Hamiltonian Hn . 
Its generalization for a k -dependent SO coupling con- 
stant, aji = ctii{k'^) , should split the O -function into 
two terms. 

It is a remarkable fact that the characteristic length 
ka^ found here for the spin production coincides with 
the Dyakonov-Perel spin-relaxation length. Indeed, the 
spin relaxation time is [akpYrp/h^ , and there- 

fore the spin diffusion length is Ls ~ vp{TpTs) ~ A:~^ . 
This coincidence indicates strong entanglement of the 
production and relaxation processes in the inhomoge- 
neous spin transport. Moreover, when q 2ka and the 
low-frequency wing of the interbranch absorption softens, 
(lj — W1/2) , the inter- and intrabranch absorption 
merge. 



Discussion and Conclusions 

The complexity of the spin transport in media with 
spin-orbit interaction produced controversies in the re- 
sults of different theoretical groups that have not been re- 
solved yet |23|2l|23,'28','2^,HlMH|5S|5l|5l|5l|5i. 

The unsettled problems include some analytical results, 
their physical interpretation, and the comparison of the 
analytical and numerical results. These problems seem to 
stem from the competition of different terms that is sim- 
ilar to the cancelations discussed in Sec. III.B. However, 
the proper treatment of such a competition for nonequi- 
librium systems is much more demanding. 
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In this paper, relation between spin currents and the 
SO contribution to the spontaneous electric polarization, 
as well as their connection to the dielectric function of 
noncentrosymmetric crystals, have been established. The 
effect of the crystal symmetry on the background (equi- 
librium) spin currents has been investigated, and the im- 
plications for spin currents driven by external electric 
fields have been discussed. Spin-response to an inhomo- 
geneous time-dependent electric field has been found. 

As distinct from spin currents that have no rigor- 
ous theoretical justification and manifest itself in di- 
electric polarization, spin magnetization is a well de- 
fined quantity. From this standpoint, spin magnetiza- 
tion and electric polarization (electric current) can be 
considered as basic variables. They are coupled through 
spin-orbit interaction and their propagation is a collec- 
tive effect [2^ • Such an approach has the advantage 
of operating in terms of measurables, i.e., the quantities 
that can be controlled experimentally, and well matches 
Maxwellian electrodynamics. 

In the practical aspect, an electric field E can produce 
spin magnetization only by breaking the time-inversion 
symmetry. When B = , it is possible either due to dis- 
sipation or by employing transients. The latter approach 
seem to have the advantage of lesser losses. Similarly, 
small structure sizes or large gradients are needed for a 
low-losses ballistic spin transport. Therefore, the geome- 
tries including transients and sharp gradients facilitate 
spin magnetization and ballistic spin fluxes. 
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